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Abstract — The Elmore delay is fast becoming ineffective for 
deep submicron technologies, and reduced order transfer 
function delays are impractical for use as early-phase design 
metrics or as design optimization cost functions. This paper 
describes an accurate approach for fitting moments of the 
impulse response to probability density functions so that delay 
and slew metric can be estimated accurately at an early 
physical design stage. PERI (Probability distribution function 
Extension for Ramp Inputs) technique has been used that 
extends the delay and slew metrics for step inputs to the more 
general and realistic non-step or ramp inputs. The accuracy 
of the proposed model is justified by the results obtained from 
the proposed model and that of SPICE simulations. 

Index Terms — Moment Matching, Delay Calculation, Slew 
Calculation, Inverse Gamma Distribution, VLSI. 

I. Introduction 

The advent of sub-quarter-micron IC technologies has 
forced dramatic changes in the design and manufacturing 
methodologies for integrated circuits and systems. The 
paradigm shift for interconnect which was once considered 
just a parasitic but can now be the dominant factor to determine 
the integrated circuit performances. It results the greatest 
impetus for the change in existing methodologies. As 
integrated circuit feature sizes continue to scale well below 
0. 1 8 microns, active device counts are reaching hundreds of 
millions [3]. This paper proposes an extension of Elmore's 
approximation [1] to include matching of higher order moments 
of the probability density function. Specifically, using a time- 
shifted incomplete gamma function approximation [2] for the 
impulse responses of RC trees, the three parameters of this 
model are fitted by matching the first three central moments 
(mean, variance, skewness), which is equivalent to matching 
the first two moments of the circuit response (m , m 2 ,). 
Importantly, it is proved that such a gamma fit is guaranteed 
to be realizable and stable for the moments of an RC tree [4] . 
In this work, we used PERI technique [13] for extending the 
delay metric derived for a step input into a delay metric for a 
ramp input for RC trees that is valid over all input slews. A 
noteworthy feature of this method is that the delay metric 
reduces to the Elmore delay of the circuit under the limiting 
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case of an infinitely slow ramp, a fact first proved in [6] to 
establish the Elmore delay as an upper bound. 

The remainder of the paper is organized as follows: 
Section 2 presents the pertinent background information 
relating to the probability and circuit theory. Section 3 shows 
the proposed techniques for extending step delay and slew 
metrics to ramp inputs. Section 4 presents the simulation 
results that validate the effectiveness of the approach. Finally, 
Section 5 concludes the paper. 

II. basic Theory 

A. Moments of a Linear Circuit Response 

Let h(t) be a circuit impulse response in the time domain 
and let H(s) be the corresponding transfer function. By 
definition, H(s) is the Laplace transform of h(t) [12], 

00 

H(s)=\h(t)e- S 'dt (i) 

o 

Applying the Taylor series expansion of e ~ st about s=0 
yields, 

»M=2L^'J*'*(')* (2) 

The i th circuit-response moment, th i is defined as [12]: 

m, =^^-^t l h(t)dt (3) 
l - o 

From (2) and (3), the transfer function H(s) can be expressed 
as: 

h{s) = fh + th x s + th 2 s 2 + th 3 s 3 + ... (4) 

B. Central Moments 

It is straightforward to show that the first few central 
moments can be expressed in terms of circuit moments as 
follows [6]: 

m i 2 ^ _m,m, „m. 3 

M<,=m ,M 1 =0,M2= 2m 2 . M, = -°m 3 + 6 -2— (5) 

m a m a m- 

Unlike the moments of the impulse response, the central 
moments have the geometrical interpretations: 
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/u Q is the area under the curve. It is generally unity, or else a 

simple scaling factor is applied. ju 2 is the variance of the dis- 
tribution which measures the spread or the dispersion of the 
curve from the center. A larger variance reflects a larger spread 

of the curve. ^ is a measure of the skewness of the distribu- 
tion. 

III. Proposed Work 

A. Calculation of the Delay Metric 

Elmore's distribution interpretation can be extended 
beyond simply estimating the median by the mean if higher 
order moments can be used to characterize a representative 
distribution function. Once characterized, the delay can be 
approximated via closed form expression or table lookup of 
the median value for the representative distribution family. 
The inverse gamma distribution is a two parameter continuous 
distribution [4] . The PDF of Inverse Gamma's distribution [4] 
is depicted in Figure 1. The inverse gamma distribution's 
probability density function is defined over x > 0, 

Where a isthe shape parameter andp isthe scale parameter. 

Mean for the inverse gamma distribution is given by, 

Mean{E[x]) = (13) 
Variance of the inverse gamma distribution is given by, 

Variance = @- (14) 

(a -if (a -2) K ' 

Mode of the inverse gamma distribution is given by, 

Mode = -P- (15) 
a + 1 

The parameters a and P of inverse gamma distribution can be 
presented in terms of moments. Moment generating function 
of inverse gamma distribution is given as, 

When, n=l; e[x] = -^ (17) 

When,n=2; (g _^ _ 2) (18) 

Hence we can write, 

OTl =-JL (19) 

m ^(dk^) (20) 

From (19) and (20), 

(a 2 -3a + l}n 2 = (a-l) 2 m, 2 (21) 



Prabillliry dinnly funni™ 




Figure 1. Inverse Gamma Distribution PDF. 

After solving above equation, 

a 2 {in 2 - >n\)+ a^lm 2 — 3m 2 )+(m 2 — m 2 )=0 
By solving (22), we get, 



- (2m 1 — 3m-, )± -J (2m, 2 - 3m-, ) - 4(m 2 - tn} \2rn-, — m} ) 

a= — A A 

4\m 2 — m l J 

-[2m? -3/tit )±J\4mf +9iUj -\2m}m 2 - 8m| + Bmfnu +4m?m 2 -4m, 4 ) 

a=^ -J— ^ i — '- 

4lm 2 — m x I 

By taking positive sign, we get for a=l/2, 

4m, — 2m, 2 



(22) 



(23) 



(24) 



4{m 2 —m\ j 

From (19) and (24), we can write for P as, 

4(m 2 — m 2 j 

The median of the inverse gamma distribution is defined as 
[9], 



(25) 



Median = — \Mode + 2Mean] 
Substituting (13) and (15), we get, 



Median - 



2/3 , P 
a -I a + 1 



After solving, 

(3a +l)p 



Median - 



(26) 



(27) 



(28) 



3(a + lX«-l) 

Now substituting the values of a and P from (24) and (25) in 
(28), we have, 



Median - 



4lm 2 — m 1 
2\m-, - m. 



2m~\ 



{ + 1 



4(m^ - mj 2 ] 
21m 7 - m. 



4lm 2 — m, j j\ 4(m 2 - mf I 

After simplification of the above equation, we get, 
m 1 (8yw 2 -5m! 2 ) 



50% Delay = Median ■ 



(29) 



(30) 
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3(4m 2 -3m, 2 ) 

The above equation (30) is the proposed closed form delay 
expression for generalized RC global interconnects using 
inverse gamma distribution based on moment generating 
technique. 

B. Calculation for the Slew Metric 

The step response of an RC circuit is a cumulative density 
function (CDF) [7]. The RC response is considered as a single- 
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pole exponential waveform and can be modeled as, 

t 

t > , If h(t) satisfies the following conditions: 



= 1- 



< h(t) < 1 and 

Now, let T 



/(f) = i 



(31) 



him him 
f(t) = 0, 
t -> -oo t - 

and T be 10% and 90% delay points, 



respectively. Matching to these points to the CDF yields, 



0.1 = l-e 



0.9 = l-e 



From (32) and (33) we have, 



(32) 
(33) 

(34) 
(35) 

Using (34) and (35), we define the inverse gamma slew metric 
(we call this metric as IGSM) as, 



/,,, - : lnl — | = 0.1053/? 



IGSM=T HI - T LO =2.1976/? 



(36) 



Using equation (24) and (25), we can write the closed form 
expression of the slew metric in terms of first two circuit 
moments as 



IGSM 



4.3952m? 



(37) 



From the above derived equation (37) for the slew metric for 
the on-chip interconnect using inverse gamma distribution 
function. 

C. Delay Metric for Ramp Input 

Let u(s) = -mj is the Elmore delay and M(S) is the step 
delay metric as given (30). The delay estimation for the ramp 
response [13] is given by, 

D(R)=(l-a)u(s)+aM(S) (38) 
Where a denotes the constant and is given as, 



2m2 - 1 
2m 2 - mf 



T 

12 ) 



(39) 



Where 0<T<" is the slope of the ramp input. 
From (30) and (38), we get, 

D(R) = -m x +a\ --- 

1.3m, 

Substituting the value of a from (39) in (40) we get, 



D(R) = -m, 



2m 2 
. 3m, 



2m 2 — /"i 
2m 2 — mf 



T 

"12 ; 



(40) 



(41) 



Equation (41) is the delay metric using inverse gamma 
distribution function for ramp input. 

D. Slew Metric for Ramp Input 

The output slew is the root-mean square of the step slew 
and input slew [ 14] . For ramp slew 

Slew(R) = ^Slew 2 (S) + Slew 2 (I) (42) 
Further, (42) exhibits the right limiting behavior: as 
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Slew(I)^x, , we have Slew(R) -> &> and as Slew(I) -> , we 
have Slew(R) Slew(S) . Where Slew(S) is the step slew metric 
which is given by (37) and Slew (I) is the input slew which is 
given as, 

Slew 2 (I)=T 2 /12 (43) 
From (37), (42) and (43), we get 



Slew(R) ■■ 



4.3952m 
4^m 2 —m { 



12 



(44) 



The above equation (44) is the slew metric equation for the 
inverse gamma distribution function for ramp input. 

IV. Simulation Results and Discussions 

We have implemented the proposed delay estimation 
method using inverse gamma distribution and applied it to 
widely used actual interconnect RC networks as shown in 
Figure-2. For each RC network source we put a driver, where 
the driver is a voltage source followed by a resister. We classify 
the 2224 sinks as it was taken in PERI [14] 

• 1 187 far-end sinks have delay greater or equal to 75% of the 
maximum delay to the furthest sink in the net. 

• 670 mid-end sinks which have delay between 25% and 75% 
of the maximum delay and, 

• 367 near-end sinks which have delay less than or equal to 
25% of the maximum delay. 




Figure 2. An RC Tree Example 

We compare the obtained average, minimum, maximum values 
along with standard deviation for delay from PERI with those 
found using our proposed model. The comparative results 
are summarized in Table 1 . 

Table I. 

Comparison of the Proposed Delay Model with Peri 
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Table I. 

Comparison of the Proposed Slew Model with Peri 
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We compare the obtained average, minimum, maximum val- 
ues along with standard deviation for slew from PERI with 
those found using our proposed model. The comparative 
results are summarized in Table 2. For the final step for experi- 
ments, we have estimated the effect of input slew on the 
delay and slew estimation for the same seven node RC net- 
work as shown in the figure 3. we have chosen 10/90 input 
slew values of 200 and 500 picoseconds and delays and slews 
were estimated for each node. The comparative result of our 
proposed model with RICE [14] are shown in the Table 3 and 
Table4 for delay and slew, respectively. 



Table III. Comparison of the Proposed Slew Model with Peri 
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Table I. 

Comparison of the Proposed Slew Model with Peri 
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V. Conclusions 

In this paper, we have extended the method to estimate 
the interconnect delay metric for high speed VLSI designs 
from step input to ramp input. We have used inverse gamma 
probability distribution function to derive our metric. Our 
model has Elmore delay as upper bound but with significantly 
less error. The novelty of our approach is justified by the 
comparison made with that of the results obtained by SPICE 
simulations. 
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